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Abstract. The cooling of compact isolated objects for different values of the gravitational 
mass has been simulated for two alternative assumptions. One is that the interior of the 
D ' star is purely hadronicp] and second that the star can have a rather large quark core It 

has been shown that within a nonlocal chiral quark model the critical density for a phase 
00 i transition to color superconducting quark matter under neutron star conditions can be low 

{SJ ' enough for these phases to occur in compact star configurations with masses below 1.3 Mq. 

For a realistic choice of parameters the equation of state (EoS) allows for 2SC quark matter 
with a large quark gap ~ 100 MeV for u and d quarks of two colors that coexists with 
normal quark matter within a mixed phase in the hybrid star interior. We argue that, if 
in the hadronic phase the neutron pairing gap in 3P2 channel is larger than few keV and 
the phases with unpaired quarks are allowed, the corresponding hybrid stars would cool too 
fast. 

Even in the case of the essentially suppressed 3P2 neutron gap if free quarks occur for 
M < 1.3 Mq, as it follows from our EoS, one could not appropriately describe the neutron 
\ star cooling data existing by today. 

It is suggested to discuss a "2SC+X" phase, as a possibility to have all quarks paired in 
two-flavor quark matter under neutron star constraints, where the X-gap is of the order of 
10 keV - 1 MeV. Density independent gaps do not allow to fit the cooling data. Only the 
presence of an X-gap that decreases with increase of the density could allow to appropriately 
fit the data in a similar compact star mass interval to that following from a purely hadronic 
c/2 . model. 

> 

X : 1. INTRODUCTION 

The "standard" scenario of neutron star cooling is based on the main process 
responsible for the cooling, which is the modified Urea process (MU) nn — > npeu 
calculated using the free one pion exchange between nucleons, see (3. However, 
this scenario explains only the group of slow cooling data. To explain a group of 
rapid cooling data "standard" scenario was supplemented by one of the so called 
"exotic" processes either with pion condensate, or with kaon condensate, or with 
hyperons, or involving the direct Urea (DU) reactions, see 0, HI and refs therein. 
All these processes may occur only for the density higher than a critical density, 
(2-h6) Uq, depending on the model, where n is the nuclear saturation density. An 
other alternative to "exotic" processes is the DU process on quarks related to the 
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phase transition to quark matter. 

Recently, the cooling of neutron stars has been reinvestigated within a purely 
hadronic model i.e., when one suppresses the possibility of quark cores in 
neutron star interiors. We have demonstrated that the neutron star cooling data 
available by today can be well explained within the "nuclear medium cooling" 
scenario, cf. 0,0], i.e., if one includes medium effects into consideration. In the 
"standard plus exotics" scenario for hadronic models the in-medium effects have 
not been incorporated, see [1,0,ED|- Recently called this approach the "minimal 
cooling" paradigm. Some papers included an extra possibility of internal heating 
that results in a slowing down of the cooling of old pulsars, see and Refs. 
therein. 

The necessity to include in-medium effects into the neutron star cooling is based 
on the whole experience of condensed matter physics, see 0,0, El- The relevance 
of in-medium effects for the neutron star cooling problem has been shown by 
0, EH EH EE E3 an d the efficiency of the developed "nuclear medium cooling" 
scenario for the description of the neutron star cooling was demonstrated within 
the cooling code by :6<] and then by 

Each scenario puts some constraints on dense matter equation of state (EoS). In 
particular the density dependencies of the asymmetry energy and the pairing gaps 
are the regulators of the heat production and transport. The former dependence is 
an important issue for the analysis of heavy ion collisions especially within the new 
CBM (compressed baryon matter) program to be realized at the future accelerator 
facility FAIR at GSI Darmstadt. 

The density dependence of the asymmetry energy also determines the proton 
fraction in neutron star matter and thus governs the onset of the very efficient 
direct Urea (DU) process. The DU process, once occurring, would lead to a very 
fast cooling of neutron stars. Within the "standard + DU" scenario the transition 
from slow cooling to the rapid cooling occurs namely due to the switching on the 
DU process. Thus the stars with M < cool down slowly whereas the stars with 
the mass only slightly above cool down very fast. Since it is doubtful that 

many neutron stars belonging to an intermediate cooling group have very similar 
masses, from our point of view such a scenario seems unrealistic, cf. The 
modern EoS of the Urbana-Argonne group (HI allows for the DU process only for 
very high density n > 5n (where n is the saturation nuclear density) that relates 
to the neutron star masses M > ~ 2 M & . Thus, using mentioned Urbana- 
Argonne based EoS and the "standard +DU" scenario one should assume that the 
majority of experimentally measured cooling points relates to very massive neutron 
stars that seems us still more unrealistic. 

The assumption about the mass distribution can be developed into a more 
quantitative test of cooling scenarios when these are combined with population 
synthesis models. The latter allow to obtain Log N - Log S distributions for nearby 
coolers which can be tested with data from the ROSAT catalogue |2(|. Analysis 
j2(^ has supported ideas put forward in 

At high star masses the central baryon density exceeds rather large values 
n > 5n . At these densities exotic states of matter as, e.g., hyperonic matter or 
quark matter perhaps are permitted. Ref. [H( argued that the presence of the 



quark matter in massive compact star cores is a most reliable hypothesis. 

The possibility of the existence of neutron stars with large quark matter cores is 



also not excluded 0, H2, IE 24 1 . In the quark matter the DU process yielding the 



rapid cooling may arise on interacting but unpaired quarks [25j 

In this review we want to sketch a scenario for the cooling of hybrid stars. 



2. STRUCTURE OF HYBRID NEUTRON STARS 

In describing the hadronic part of the hybrid star, as in ij , we exploit a modification 
of the Urbana-Argonne V18 + 5v + UIX* model of the EoS given in [lflj . which 
is based on the most recent models for the nucleon-nucleon interaction with the 
inclusion of a parameterized three-body force and relativistic boost corrections. 
Actually we continue to adopt an analytic parameterization of this model by 
Heiselberg and Hjorth- Jensen [26], hereafter HHJ. 

The HHJ EoS fits the symmetry energy to the original Argonne V18 + 8v + UIX* 
model in the mentioned density interval yielding the threshold density for the DU 
process n^ u - 5.19 n (M C DU ~ 1.839 M ). 

We employ the EoS of a nonlocal chiral quark model developed in for the case 
of neutron star constraints with a 2-flavor color superconductivity (2SC) phase. It 
has been shown in that work that the Gaussian formfactor ansatz leads to an early 
onset of the deconfinement transition and such a model is therefore suitable to 
discuss hybrid stars with large quark matter cores p8^ . 

The quark-quark interaction in the color anti-triplet channel is attractive driving 
the pairing with a large zero-temperature pairi ng g ap A ~ 100 MeV for the quark 
chemical potential fi q ~ (300^-500) MeV, cf. |2flL l3oT|. for a review see [Hj[ and 
references therein. The attraction comes either from the one-gluon exchange, or 
from a non-perturbative 4-point interaction motivated by instantons [I]]], or from 
non-perturbative gluon propagators |3^ |. 

There may also exist a color-flavor locked (CFL) phase for not too large 
values of the dynamical strange quark mass or large values of the baryon chemical 
potential In this phase the all quarks are paired. However, the 2SC phase 
occurs at lower baryon densities than the CFL phase, see [sHE^l- For applications 
to compact stars the omission of the strange quark flavor is justified by the fact 
that chemical potentials in central parts of the stars do barely reach the threshold 
value at which the mass gap for strange quarks breaks down and they appear in 
the system 

Following Refs |38| we omit the possibility of the hadron-quark mixed phase 
and found a tiny density jump on the phase boundary from n)? adr ~ 0.44 fm~ 3 to 
n quMk~0.46 fm' 3 . 

In Fig. we present the mass-radius relation for hybrid stars with HHJ EoS vs. 
Gaussian nonlocal chiral quark separable model (SM) EoS. Configurations for SM 
model, given by the solid line, are stable, whereas without color super conductivity 
("HHJ-SM without 2SC") no stable hybrid star configuration is possible. In the 
case "HHJ-SM with 2SC" the maximum neutron star mass proves to be 1.793 M & . 
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FIGURE 1. Mass - radius relations for compact star configurations with different EoS: purely 
hadronic star with HHJ EoS (dashed line), stable hybrid stars with HHJ - Gaussian nonlocal 
chiral quark separable model (SM) with 2SC phase (solid line) and with HHJ - SM, without 2SC 
phase (dash-dotted line) . Data for two sources are also indicated, see [3^, |4jj . 

Additionally, within the "HHJ-SM with 2SC" phase we will allow for the pos- 
sibility of a weak pairing channel for all the quarks which were unpaired, with 
typical gaps A x ~ 10 keV -r-1 MeV, as in the case of the CSL pairing channel, see 
fill H2]. Since we don't know yet the exact pairing pattern for this case, we call 
this hypothetical phase "2SC+X" . In such a way all the quarks get paired, some 
strongly in the 2SC channel and some weakly in the X channel. 

3. COOLING 

We compute the neutron star thermal evolution adopting our fully general rela- 
tivists evolutionar y co de. This code was originally constructed for the description 
of hybrid stars by |2jJ . The main cooling regulators are the thermal conductivity, 
the heat capacity and the emissivity. In order to better compare our results with 
results of other groups we try to be as close as possible to their inputs for the 
quantities which we did not calculate ourselves. Then we add inevitable changes, 
improving EoS. 

The density n ~ 0.5 -j- 0.7 n is the boundary of the neutron star interior and the 
inner crust. The latter is constructed of a pasta phase discussed by see also 
recent works of pHH^. 

Further on we need the relation between the crust and the surface temperature 
for neutron star. The sharp change of the temperature occurs in the envelope. 



3.1. Cooling Evolution of Hadronic Stars 



Here we will shortly summarize the results on hadronic cooling. 

In framework of "minimal cooling" scenario the pair breaking and formation 
(PBF) processes may allow to cover an "intermediate cooling" group of data (even 
if one artificially suppressed medium effects) [6|. These processes are very efficient 
for large pairing gaps, for temperatures being not much less than the value of the 
gap. 




FIGURE 2. Neutron and proton pairing gaps according to model I (thick solid, dashed and 
dotted lines) and according to model II (thin lines), see text. The ISq neutron gap is the same 
in both models, taken from |4fij . 

Gaps that we have adopted in the framework of the "nuclear medium cooling" 
scenario, see are presented in Fig. El Thick dashed lines show proton gaps 
which were used in the work of [^] performed in the framework of the "standard 
plus exotics" scenario. We will call the choice of the "3nt" model from the model 
I. Thin lines show ISq proton and 3P 2 neutron gaps from |48|, for the model AV18 
by |H| (we call it the model II). Recently has argued for a strong suppression 
of the 3P2 neutron gaps, down to values ~ 10 keV, as the consequence of the 
medium-induced spin-orbit interaction. 

These findings motivated [l| to suppress values of 3P2 gaps shown in Fig.Elby an 
extra factor f(3P 2 ,n) = 0.1. Further possible suppression of the 3P 2 gap is almost 
not reflected on the behavior of the cooling curves. 

Contrary to expectations of Ji?} a more recent work of [E(| argued that the 
3P2 neutron pairing gap should be dramatically enhanced, as the consequence of 
the strong softening of the pion propagator. According to their estimate, the 3P2 
neutron pairing gap is as large as 1 -j- 10 MeV in a broad re gion of densities, see 
Fig. 1 of their work. Thus results of calculations of [47| and |50j ] , which both had 
the same aim to include medium effects in the evaluation of the 3P2 neutron gaps, 
are in a deep discrepancy with each other. 

• Including superfluid gaps we see, in agreement with recent microscopic findings 
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FIGURE 3. Fig. 21 of Gaps are from Fig.|2]for model II. The original 3P2 neutron pairing 
gap is additionally suppressed by a factor /(3iVn) = 0.1. The T s — T m relation is given by "our 
fit" curve of Fig. 4 in pj . For more details see [lj . 

of that 3P 2 neutron gap should be as small as 10 keV or less. So the 
"nuclear medium cooling" scenario of [l[ supports results of j4?J and fails to 
appropriately fit the neutron star cooling data at the assumption of a strong 
enhancement of the 3P 2 neutron gaps as suggested by [Ell ]. 

• Medium effects associated with the pion softening are called for by the data. 
As the result of the pion softening the pion condensation may occur for 
n > (n > 3uq in our model). Its appearance at such rather high densities 
does not contradict to the cooling data (see Fig. |3J), but also the data are 
well described using the pion softening but without assumption on the pion 
condensation. This also means that the DU threshold density can't be too 
low that puts restrictions on the density dependence of the symmetry energy. 
Both statements might be important in the discussion of the heavy ion collision 
experiments. 

• We demonstrated a regular mass dependence: for the neutron star masses 
M > 1 M & less massive neutron stars cool down slower, more massive neutron 
stars cool faster. 



3.2. Cooling Evolution of Hybrid Stars with 2SC Quark 

Matter Core 

For the calculation of the cooling of the quark core in the hybrid star we use 
the model j24|]. We incorporate the most efficient processes: the quark direct Urea 
(QDU) processes on unpaired quarks, the quark modified Urea (QMU), the quark 
bremsstrahlung (QB), the electron bremsstrahlung (EB), and the massive gluon- 
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FIGURE 4. Model I. Cooling curves for hybrid star configurations with Gaussian quark matter 
core in the 2SC phase. The labels correspond to the gravitational masses of the configurations in 
units of the solar mass. 

photon decay (see [22| ) • Following [H^ we include the emissivity of the quark pair 
formation and breaking (QPFB) processes. The specific heat incorporates the quark 
contribution, the electron contribution and the massless and massive gluon-photon 
contributions. The heat conductivity contains quark, electron and gluon terms. 

The calculations are based on the hadronic cooling scenario presented in Fig. |H] 
and we add the contribution of the quark core. For the Gaussian form-factor the 
quark core occurs already for M > 1.214 M according to the model (22J, see Fig. 
^ Most of the relevant neutron star configurations (see Fig. |3J) are then affected 
by the presence of the quark core. 

First we check the possibility of the 2SC+ normal quark phases Fig. 0] 

The variation of the gaps for the strong pairing of quarks within the 2SC phase 
and the gluon-photon mass in the interval A,m 9 _ 7 ~ 20 -=-200 MeV only slightly 
affects the results. The main cooling process is the QDU process on normal quarks. 
We see that the presence of normal quarks entails too fast cooling. The data could 
be explained only if all the masses lie in a very narrow interval (1.21 < M/M Q < 1.22 
in our case). In case of the other two crust models the resulting picture is similar. 

The existence of only a very narrow mass interval in which the data can be 
fitted seems us unrealistic as by itself as from the point of view of the observation 
of the neutron stars in binary systems with different masses, e.g., M B i9i3+i6 — 
1.4408 ±0.0003 M and M J07 37-3039B ^ 1.250 ±0.005 M , cf. @. Thus the data 
can't be satisfactorily explained. 

We first check the case Ax to be constant. For the ~ 1 MeV cooing is too 
slow 0. It is true for all three crust models. Thus the gaps for formerly unpaired 
quarks should be still smaller in order to obtain a satisfactory description of the 
cooling data. 



For the Ax = 30 keV the cooling data can be fitted but have a very fragile 
dependence on the gravitational mass of the configuration. Namely, we see that 
all data points, except the Vela, CTA 1 and Geminga, correspond to hybrid stars 
with masses in the narrow interval M = 1.21 -r 1.22 M & 

Therefore we would like to explore whether a density-dependent X-gap could 
allow a description of the cooling data within a larger interval of compact star 
masses. 

We employ the ansatz: X-gap as a decreasing function of the chemical potential 

A x (/i) = A c exp[-a(/i-/i c )//i c ] , (1) 

where the parameters are chosen such that at the critical quark chemical potential 
fi c = 330 MeV for the onset of the deconfinement phase transition the X-gap has its 
maximal value of A c = 1.0 MeV and at the highest attainable chemical potential 
AWx = 507 MeV, i.e. in the center of the maximum mass hybrid star configuration 
it falls to a value of the order of 10 keV. We choose the value a = 10 for which 
Ax(/i max ) = 4.6 keV. In Fig. El we show the resulting cooling curves for the gap 
model II with gap anzatz^ which we consider as the most realistic one. 

We observe that the mass interval for compact stars which obey the cooling 
data constraint ranges now from M = 1.32 M Q for slow coolers up to M = 
1.75 M & for fast coolers such as Vela, cf. with that we have found with the purely 
hadronic model |]J with different parameter choices. Note that according to a 
recently suggested independent test of cooling models [2(j by comparing results of 
a corresponding population synthesis model with the Log N - Log S distribution of 
nearby isolated X-ray sources the cooling model I did not pass the test. Thereby 
it would be interesting to see whether our quark model within the gap ansatz II 
could pass the Log N - Log S test. 
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FIGURE 5. Cooling curves for hybrid star configurations with Gaussian quark matter core 
in the 2SC phase with a density dependent pairing gap according to Eq. JQ) for model II. 



4. CONCLUSION 



• Within a nonlocal, chiral quark model the critical densities for a phase transi- 
tion to color superconducting quark matter can be low enough for these phases 
to occur in compact star configurations with masses below 1.3 M & . 

• For the choice of the Gaussian form-factor the 2SC quark matter phase arises 
at M ~ 1.21 M Q . 

• Without a residual pairing the 2SC quark matter phase could describe the 
cooling data only if compact stars had masses in a very narrow band around 
the critical mass for which the quark core can occur. 

• Under assumption that formally unpaired quarks can be paired with small 
gaps Ax < 1 MeV (2SC+X pairing), which values we varied in wide limits, 
only for density dependent gaps the cooling data can be appropriately fitted. 

So the present day cooling data could be still explained by hybrid stars, however, 
when assuming a complex pairing pattern, where quarks are partly strongly paired 
within the 2SC channel, and partly weakly paired with gaps Ax < 1 MeV, being 
rapidly decreasing with the increase of the density. 

It remains to be investigated which microscopic pairing pattern could fulfill the 
constraints obtained in this work. Another indirect check of the model could be 
the Log N - Log S test. 
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